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ABSTRACT 

High-resolution optical spectra of the ultraluminous X-ray source NGC 5408 X-1 show a broad 
component with a width of ~750 km/s in the Hen and H/3 lines in addition to the narrow component 
observed in these lines and [O ill] . Reanalysis of moderate- resolution spectra shows a similar broad 
component in the Hen line. The broad component likely originates in the ULX system itself, probably 
in the accretion disk. The central wavelength of the broad Hell line is shifted by 252 ± 47 km/s 
between the two observations. If this shift represents motion of the compact object, then its mass is 
less than -1800 Mq. 

Subject headings: black hole physics — galaxies: individual(NGC 5408) — X-rays: binaries 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are variable off- 
nuclear X-ray sources with luminosities exceeding the 
Eddington luminosity of a 20 Mo compact object, as- 
suming isotropic e mission ([Colbert fc Mushotzkv I [l999l : 
iKaaret et al.ll200l|) . Irregular variability, on time scales 
from seconds to years, suggests that ULXs contain ac- 
creting compact objects. Intermediate mass black holes 
would be required to produce the inferred luminosities, 
but ULXs may, instead, accrete at super Eddington rates 
or be beamed, mechanically or relativistically. 

NGC 5408 X-1 is one of the best intermediate mass 
black hole candidates because it powers a radio nebula 
requi r ing an extremely energetic outflow (jKaaret et al.l 
120031 iSoria et al.l 120061 : iLang et al.l l2007| ) and a pho- 
toionized nebula req uiring an X -ray luminosity above 
3 X lO^^ergs"^ (Kaa iet fc Corb el 2009). Also, quasi- 
periodic X-ray oscillations at low frequencies suggest a 
high compact object mass (jStrohmaver et al. H2007). 

The optical coun t erpart to NGC 5408 X-1 was identi- 
fied by iLang et all (I2007D and op tical spectra were ob- 
tained bv K aaret fc Corbel (|2009[ ). The optical spectra 
had no absorption lines suggesting the emission is not 
dominated by the companion star. The observed contin- 
uum emission may arise from a nebula or reprocessing 
of X-rays in an accretion disk. The optical spectrum is 
dominated by emission lines, including forbidden lines 
which must be produced in a low density environment 
such as a nebula. Several high excitation lines were de- 
t ected indicating that the nebula is X-ray photoionizcd. 

IKaaret fc Corbell (|2009[ ) found that the Heii line from 
NGC 5408 X-1 was broader than the forbidden lines. 
Permitted lines produced in the high-density environ- 
ment of an accretion disk can be broad, reflecting the 
distribution of velocities within the optical emitting re- 
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gions of the disk. Furthermore, since the accretion disk 
moves with the compact object, the line velocity shifts 
may provide a means to constrain the compact object 
mass (jHutchings et al.lll987HSoria et al.|[T998il . 

To study the Hen line profile of NGC 5408 X-1 in 
more detail, we obtained new observations using the 
FORS-2 spectrograph on the European Southern Obser- 
vatory Very Large Telescope (VLT) with a high resolu- 
tion gris m and reanalyzed our previous FORS-1 obser- 
vations (jKaaret fc Corbell I2009D . The observations and 
data reduction are described in §2. The results are pre- 
sented in §3 and discussed in §4. 

2. OBSERVATIONS AND ANALYSIS 

FORS-2 observations of NGC 5408 X-1 were obtained 
on 12 April 2010 using the GRIS_1200B and GRIS_1200R 
grisms with a slit width of 1.0" covering the spectral 
range 3660-5110 A and 5750-7310 A with dispersion 
0.36 A pixel^^ and 0.38 A pixel"^ and spectral resolu- 
tion A/AA = 1420 and A/AA = 2140 at the central wave- 
length, respectively. The observation block (OB) con- 
sisted of three 849 s exposures with a 12 pixel offset along 
the spatial axis between successive exposures. CCD pix- 
els were binned for readout by 2 in both the spatial 
and spectral dimensions. We also reanalyzed all six OBs 
of th e previous FORS-1 observations (Kaarct ~fc Corbell 
l2009f l. hereafter the low resolution data (LRD), taken 
using the GRIS_600B grism which has a spectral resolu- 
tion of A/AA = 780 at the central wavelength and with 
three shifted exposures per OB. The average seeing for 
our new observations was 0.72 and 0.62 arcsecond for the 
blue and red spectra, respectively. The average seeing of 
the six OBs of the LRD were 0.87, 0.82, 0.96, 1.28, 0.64, 
0.57 arcsecond, respectively. 

Data reduction was carried out usi ng the Image Reduc- 
tion and Analysis Facility (IRAPfl (jTodv II1993I ). First, 
we created bias and flat-field images, then applied these 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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to correct the spectrum images. The three exposures in 
each OB were ahgned then averaged to ehminate bad 
pixels and cosmic rays using the imcombine task with 
the ccdclip rejection algorithm. 

As the continuum emission of the ULX counterpart 
is faint, we could not tr ace its spectrum. Follow- 
ing |KaareF^£o^e|j200^}j_wre used the bright nearby 
star at 2MASS (jSkruitske et al.|[2006l ) position q;j2ooo = 
14''03"'18.^97,(5j2ooo = -41°22'56."6 as a reference 
trace. The trace position on the spatial axis varied less 
than half a pixel along the whole length of the disper- 
sion axis. The trace for the ULX counterpart was cen- 
tered on the Hen A4686 emission line profile. The small- 
est possible trace width, 2 pixels corresponding to 0.5", 
was used to best isolate the ULX emission from the neb- 
ular emission. Background subtraction was done with 
a trace close by. The HgCdHeNeA lamp and standard 
star LTT7379 were used for wavelength and flux calibra- 
tion. An atmospheric extinction correction was applied 
using the IRAF built-in Cerro Tololo Inter- American Ob- 
servatory (CTIO) extinction tables. To estimate the 
reddening, we used the Balmer decrement of H5/H/3, 
we find E(B - V) ^ 0.08 ± 0.03 in agreement with 
iKaaret fc Corbell ()2009[ ). We corrected for reddenin g us- 
ing the extinction curve from lCardelh et all (|1989D with 
Rv = 3.1. 

To study the kinematics, we need to characterize the 
instrumental resolution in order to obtain intrinsic line 
widths. After applying the dispersion correction to the 
lamp spectrum, we measured the full width at half maxi- 
mum (FWHM) of several lines, excluding saturated ones, 
by fitting Gaussians with the IRAF splot subroutine. 
The instrumental FWHM was 2.24 A and 5.08 A for the 
high and low resolution data, respectively. The error on 
the instrumental FHWM was estimated by finding the 
standard deviation of the FWHM for several different 
lines. 

For the Hen, H/3 and [Oni] emission lines, see 
Fig. 1. and Table [H we first fitted the continuum 
with a second order polynomial to a region around 
each line excluding the line itself by visual examina- 
tion. We estimated the measurement errors by calcu- 
lating the root mean square deviation of the data in 
the same region. Then, we performed a non-linear least 
squares fit using the LMFIT subroutine of the Interactive 
Data Language version 7.0 and based on "MRQMIN" 
(jTeukolskv. Vettering fc Flannerv |[l992l ). We fitted the 
line profiles iteratively, first using one Gaussian which 
converged on the narrow component, then using a sum 
of two Gaussians with initial parameters adjusted to 
achieve convergence. All six parameters in the two Gaus- 
sian fit were free to vary. The errors on the parame- 
ters were calculated by the fitting routine in a way that 
the uncertainty for the ith parameter derives from the 
square-root of the corresponding diagonal element of the 
covariance matrix. The intrinsic line width was calcu- 
lated assuming the measured line width is the quadra- 
ture sum of the intrinsic and instrumental widths and 
the error on the intrinsic line width included a term for 
the uncertainty in the instrumental FHWM. 

For the Ha line, we fitted the sum of four Gaussians, 
because the [Nn] lines lie on the red and blue parts of 
the line wing. Initial fits to the Ha and red [Nn] lines 
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Fig. 1. — Hen, H/3 and [O III] lines from the new, high resolution 
data. Results of line fits with two Gaussian components are shown 
as the separate components plotted with dashed-dot lines and the 
overall fit plotted with dashed lines (except [O ni]). 

provided initial values for a fit with four Gaussians. Be- 
cause the blue [Nn] line has very low signal to noise ra- 
tio, the widths of the two [Nn] lines were set equal, the 
wavelength offset was fixed at —35.44 A, and the am- 
plitude of the blue line was set to 1/3 of the red line 
(jOsterbrock fc Ferlandl 120061 ). 

3. RESULTS 
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Fig. 2.— Hen line from the low resolution data{OB3,OB5,OB6). 
Results of line fits with two Gaussian components are shown as the 
separate components plotted with dashed-dot lines and the overall 
fit plotted with dashed lines. 

The improved resolution of the new data clearly re- 
solves a broad component in the Heii line profile, see 
Fig. 1. and Table [1] The centroid is shifted from the 
nebular component by +0.87 ± 0.26 A in the red direc- 
tion. We also searched for broad components in other 
lines. H/3 has a broad component with a FWHM similar 
to the Heii line but shifted by —0.52 ±0.32 A towards the 
blue, rather than the red. In contrast, a single Gaussian 
provides a good fit to the forbidden [Oiii] line and there 



is no evidence for a broad component, as expected if the 
line is emitted only from the nebula. 

Then we fitted the Heii line profiles of the six OBs of 
the LRD, see Fig. 2. and Table [T] The flux variation 
of the overall line profiles correlates with the seeing, e.g. 
0B5 has the best seeing and the highest fiux. We de- 
tected a broad component in the Heii line in OBS, 0B5, 
and 0B6. We did not significantly detect a broad compo- 
nent in OBI, 0B2 and 0B4. This may be due to seeing 
or variations in the flux of the broad component. 

We note that Kaaret & Corbel (200^ reported lower 
fluxes for Heii, [Nev], and the continuum emission for 
0B4 (with by far the worst seeing) as compared to the 
other OBs, while the other line fluxes remained relatively 
constant. Our new analysis suggests that this is due to 
changes in the seeing. If the emitting region is smaller 
than the 0.51" slit used for the LRD, then poor seeing 
will decrease the flux through the spectrometer. If these 
emission components are enhanced close to the ULX sys- 
tem while the other line emission is uniform, then the 
poor seeing in 0B4 would produce the observed changes 
in flux. Thus, there is no evidence for temporal variabil- 
ity of the continuum or line emission. Howeve r, the sub- 
traction performed bv lKaaret fc Cor bell (|2009( ) to isolate 
continuum emission arising from near the ULX is still 
justifled; the separation of components is spatial instead 
of temporal. 

The Hen line parameters are consistent between 0B3, 
0B5, and 0B6. The wavelength shifts of the Hen broad 
component of OBS, 0B5 and 0B6 relative to the nar- 
row component are -1.82 ±0.78 A, -3.77 ± 1.44 A and 
—4.03 ± 1.70 A into the blue direction instead of the red 
as in the HRD, and are consistent within one a. We 
averaged the spectra for these three observations and fit 
the resulting line profile. The fit results are listed as 
Hen AVG in Table [TJ The shift of average line profile is 
-S.07±0.68 A. 

The Hell broad component width is consistent between 
the new and old data. The narrow component is wider 
in the old data because we do not resolve the nebular 
lines. The line fluxes are higher in the new data, most 
likely due to the wider slit. The centroids of the narrow 
component are consistent, while the wavelength shift of 
the Hen broad component between the old and new data 
is AA = S.94±0.7S A. 

Fitting the H/3 fine of the LRD, we did not get a good 
fit, due to the lack of spectral resolution and the low 
broad to narrow flux ratio. We could not fit the bluer 
Balmer lines because of their low S/N ratios. We did fit 
the Ha line in the new data. Although we do not obtain a 
good fit {xl = 4.9) because of the complicated line profile 
(i.e. the two [Nil] lines lie on the red and blue wing of 
the Ha line), we find that there is a broad component 
with a width of 19 A, while the width of the nebular 
component is 2.7 A. The [Nil] lines are narrow, with a 
typical width of 3 A, quantitatively supporting that the 
forbidden, nebular lines do not have broad components. 

4. DISCUSSION 

Our new, high-resolution spectra show narrow nebular 
lines and broad components in the Heii, H/3, and Ha 
lines. Our previous, moderate-resolution spectra show 
a broad component in the Heii line. There is no broad 
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TABLE 1 
Line fit results. 



Line 
IMF 



H/3 



Wavelength (A) 
± 0.04 



[O III] A4959 



4693.75 
4694.62 
4869.73 
4869.21 
4967.52 



Hen 0B3 
Hen OB5 
Hen OB6 
Hen AVG 



4693.90 
4692.08 
4693.80 
4690.03 
4693.97 
4689.94 
4693.89 
4690.82 



± 0.26 
± 0.01 
± 0.32 
± 0.01 
± 0.15 
± 0.77 
± 0.12 
± 1.43 
± 0.12 
± 1.70 
± 0.08 
± 0.68 



FWHM (A) 
2.54 ± 0.14 
11.98 ± 0.74 
2.28 ± 0.05 
12.49 ± 0.94 
2.48 ill 0.03 
5.15 ± 0.52 
11.09 ± 1.44 
5.81 ± 0.41 
11.75 ± 1.33 
5.56 ± 0.46 
10.82 ± 1.55 
5.50 ± 0.26 
11.47 ± 0.73 



Flux 
3.39 ± 0.13 
1.19 ± 0.10 
8.37 ± 0.11 
0.71 ± 0.07 
8.42 ill 0.08 
1.34 ± 0.19 



Velocity (km s 'y 



76.3 
752.0 

25.2 
756.9 

64.6 



3.4 
46.4 
3.3 
57.5 
1.7 



0.56 
1.74 
0.49 
1.70 
0.53 
1.58 



0.20 
0.17 
0.13 
0.23 
0.15 
0.10 



0.53 ± 0.09 



52.0 ± 33.3 
630.5 it 91.8 
179.8 it 25.7 
678.0 it 84.7 
143.7 it 29.5 
610.7 it 98.9 
135.0 it 16.6 
657.5 it 46.3 



0.94 



1.08 



1.08 



1.04 



1.03 



0.95 



1.15 



Note. — The table gives the: line identification, central wavelength, measured FWHM, flux in units of 10~^ erg cm~ s~ , intrinsic 
FWHM in velocity units, and xi of the line proflle fit. The first row for each spectral line shows the parameters of the narrow component 
and the second row shows the parameters of the broad component (except for [O m]). 



component in the [O ill] nebular lines in either the new 
or old spectra. There is still no sign of any absorption 
lines in the new spectra. 

4.1. The line emitting region 

The broad components of both Hen and H/3 have 
widths ~750 km/s, consistent with production in the ac- 
cretion disk, and are roughly Gaussian, instead of hav- 
ing P-Cygni profiles that would indicate origin in a wind. 
Following Porter (2010), we estimate the size of the line- 
emitting region, i?;e, by assuming the line-emitting gas 
is in Keplerian orbits around a compact object, thus 

Rie < GM/v^. We find Ri^ < 2.35 (jgl^) AU, which 

for a mass of 10 M0 would give an upper limit of 3.4 
Rq. This is consistent with origin of the broad Heii line 
in the accretion disk. 

The broad line components are shifted relative to the 
narrow components. In the new data, the shifts are small 
compared to the line width, -|-56 ± 17 km/s for Heii and 
—33 ± 20 km/s for H/3. These shifts are consistent only 
at the 3(7 level, which might indicate a difference in the 
spatial origin of the lines. However, this is still consis- 
tent with production of both lines within the disk since 
random motions within the disk and variation between 
the emission regions could produce shifts that are small 
compared to the line widths, as observed. 

The central wavelength of the Hen broad component 
shifts markedly between the odld and new data, AA = 
3.94 ± 0.73 A or At; = 252 ± 47 km/s. This shift is a 
substantial fraction of the line width. The shift could be 
due to random motion within the d isk, differing view- 
ing geometries ([Roberts et al.l[2010( ). or orbital motion 
of disk (and the compact object). If the shifts in the 
broad component of the Hen line are due to orbital mo- 
tion, then this would provide a means to determine the 
orbital period and would also provide a measurement of 
the mass function for the secondary star. Thus, a pro- 
gram of monitoring NGC 5408 X-1 with high- resolution 
optical spectroscopic observations will be important in 
extending our understanding of the physical nature of 
this system. 

4.2. The binary system 

In this section, we make some speculations based on 
interpretation of the shift in the broad component of 
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Fig. 3. — Properties of the binary system as a function of orbital 
period. The two solid lines show the mass of the compact object; 
the upper line is for an 120 Mq donor and the lower line is for a 
20 Mq donor (left vertical scale). The orbital separation is shown 
as dotted lines with values read on the right vertical scale; the 
upper line for an 120 M0 and the lower is for 20 donor (right 
vertical scale). The two dashed lines show the upper limit on the 
size of the line-emitting region; the upper line is for an 120 Mq 
donor and the lower line is for a 20 Mq donor (right vertical scale) . 
The intercept of the dotted and dashed lines provides a limit on 
the period and consequently on the mass of the black hole. 

the Hell line as due to orbital motion. One can express 
the mass function and the compact object mass, M^, in 
terms of the orbital period, P, the velocity excursion, 
K^, and the companion mass, M^, as 



fx 



2ttG 



Mc sin"^ i 
\ 2ttG ) 



Mr, 



(1) 



(2) 



where i is the inclination angle and G is the gravitational 
constant. From the shift of the Hen line quoted above, 
we constrain the semi-amplitude of the radial velocity 
> Av/2 = 126 ± 24 km/s. Thus, if the maximum 
mass of the companion and the orbital period are known, 
then Eq. 2 leads to an upper bound on the mass of the 
compact object. 
The binary system has a visual magnitude vq = 22.2 
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that gives an upper limit on the absolute magnitude of 
the companion of Vp = — 6.2 at a distance of 4.8 Mpc 
(|Karachentsev et al.l |2002[ ) . Unfortunately, this places 
little restriction on the companion mass as even 03V 
stars, with masses of 120 Mq, are allowed. However, 
very high mass stars are very short lived, no more than 
a few million years. There is no evidence of a dense 
stellar association near NGC 5408 X-1 and origin in the 
closest super-star cluster would require a t ransit time to 
the present location on the order of 30 Myr (^Kaar et et al.l 
11003). Thus, the companion mass is likely significantly 
lower, near 20 Mq or less similar to found fr om studies 
of the stell ar environments of other ULXs ([Grise et al.l 
[2008l[20Tll) . Fi gure 3. shows the upper bound on the 
compact object mass for donors of 120 M0 and 20 Mq 
as a function of orbital period. High black hole masses 
are excluded, except for very short periods. We note that 
the Hen line shift was the same in 0B3 versus 0B5 and 
0B6, taken one day apart, suggesting that the period is 
longer than few days. Thus, the black hole mass is likely 
below ^1800 M0. The more probable companion mass 
of 20 Mq or less would imply smaller black hole masses, 
less than 112 Mq. 

As a further constraint, we note that the orbital separa- 
tion should be larger than the size of the emitting region, 

calculated above. Assuming a circular orbit, the orbital 

1 

separation of the compact object is a = ^ ^ _ 

Figure 3. shows the the orbital separation as a function 
of period. Also shown is the size of the line-emitting 
region, Rie, versus period. Both are calculated using 
the maximum black hole mass for each period assum- 
ing an 120 Mq or 20 Mq donor. The orbital separation 
is greater than the upper limit of the size of the line- 
emitting region when the compact object mass is below 
875 M0 for an 120 Mq companion and below 128 Mq for 
a 20 Mq companion. These masses are reduced if the in- 
clination is lowered. These results suggest that the most 
probable black hole mass is at most a factor of several 
a bove the u sual stellar-mass black hole range. 

iStrohmaver 1 (|20l)9t ) proposed an orbital period of P = 
115.5 ± 4.0 days for NGC 5408 X-1, based on varia- 



tions in the X-ray emission. With P = 115.5 days, the 
mass function is fx = 24. 0± 13.4 Mq, implying a lower 
bound on the companion mass Mc > 10.6 M©. It is in- 
teresting to determine if this period is consistent with 
other constraints on the system. An orbital period of 
P = 115.5 ± 4.0 days would require a mean stellar den- 
sity of p — 1.5 X 10^^ g cm~'^ and, thus, a supergiant 
comp anion if mass transf e r proceeds via Roche-lob e over- 
flow (jStrohmaver 1 120091: iKaaret. Simet. fc Land [2006). 
In particular, late F and early G supergiants have densi- 
ties close to that required, although we caution that the 
high mass transfer rate needed to power the ULX may 
distort the spectral type of the star. Such stars have 
masses of 10-12Mq, consistent with the minimum mass 
derived from the mass function, and absolute magnitudes 
close to or below the upper limit quoted above. The stel- 
lar radii are large, up to 0.7 AU, but smaller than the 
orbital separation for this period and mass. However, a 
companion mass so close to the lower bound on the mass 
function would require a very low mass compact object. 
For a companion mass of 10-12 Mq, the compact object 
would have to be below 1 Mq, which seems unlikely. For 
a 5 Mq black hole, one would need a donor of about 
17 Mq if the system is edge on. A higher black hole 
mass or a less extreme inclination would require an even 
higher mass companion. These high companion masses 
contradict the required stellar density; there is no star 
with both Mc > 17 Mq and p ~ 1.5 x 10"^ g cm-^. 
Thus, either the orbital period is not near 115 days or 
mass transfer does not proc e ed via Roche-lobe overflow. 
We note that iFoster et al.l ()2010[ ) have suggested that 
115 day periodicity may, instead, indicate a super-orbital 
period. 
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